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Abstract 

Surface lithium abundance and rotation velocity can serve as powerful and mutually 
complementary diagnostics of interior structure of stars. So far, the processes respon- 
sible for the lithium depletion during pre-main sequence evolution are still poorly 
understood. We investigate whether a correlation exists between equivalent widths 
of Li (EW(Li)) and rotation period (P ro t) f° r Weak-line T Tauri stars (WTTSs). 
We find that rapidly rotating stars have lower EW(Li) and the fast burning of Li 
begins at the phase when star's P ro t evolves towards 3 days among 0.9M Q to 1.4M© 
WTTSs in Taurus-Auriga. Our results support the conclusion by Piau & Turch- 
Chieze about a model for lithium depletion with age of the star and by Bouvier et 
al. in relation to rotation evolution. The turn over of the curve for the correlation 
between EW(Li) and P ro t is at the phase of Zero- Age Main Sequence (ZAMS). The 
EW(Li) decreases with decreasing P roi before the star reaches the ZAMS, while 
it decreases with increasing P TO t (decreasing rotation velocity) for young low-mass 
main sequence stars. This result could be explained as an age effect of Li depletion 
and the rapid rotation does not inhibit Li destruction among low mass PMS stars. 
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1 Introduction 



T Tauri stars (TTS) are very young (< 10 8 yr), low-mass (M< 2M ), late 
spectral type (late G to M) pre-main sequence (PMS) stars, still in the process 
of gravitational contraction. They are of particular interest to stellar evolution, 
among other things because they are, if sufficiently young, in a fully convective 
stage of evolution. T Tauri stars are usually divided into two classes, namely 
Classical T Tauri stars (CTTSs) and Weak lined T Tauri stars (WTTSs). The 
former ones are associated with a circumstellar disk, from which they accrete 
material at a rate of about 10 _7 M Q yr _1 , while the later ones lack such disks. 



At the beginning of the T Tauri phase of the base of the convective zone is too 
cool to allow lithium burning. As the star evolves on the pre-main sequence 
the deep regions of convective zone temporally exceed the 7 Li burning point 
in typical stellar condition ( 2.5xl0 6 K), and surface Li can be transported to 
this region and depleted during PMS evolutionary phase. The 7 Li abundancds 
of stars therefore offer a directly insight over stellar internal structure and 
evolution as it is extremely sensitive to the appearance of the radiative core. 
The surface Li abundances and surface rotation velocity can se rve as powerf ul 



and mutually complementary diagnostics of interior structure ([Strom . 1994). 



Based on the result of Bouvier et al. ( 1993f ) from a study of the connection 



between EW(Li) and P rot of 10 CTTSs and 6 WTTSs, we concentrate on the 
connection between EW(Li) and P rot for more WTTSs. We did not consider 
CTTS, since accretion leads to a replenishment of lithium at their surface for 
CTTSs ( Bouvier et all fl993h . such that the Li-abundance cannot be used as a 



diagnostic for the interior. Also, the mass and the age for WTTSs derived from 
theoretical tracks and isochrones are more reliable than those for CTTSs. We 
performed a new and homogeneous analysis of all the Li data and P rot available 
in the early literature for WTTSs in the range of mass 0.9-1.4M Q . 



It is known that Li de pletion increases with decreasing stellar mas s , and with 
increasing metallicity (jProfhtt fc Michaudl . ll989T ). ISoderblom et all (|l99.4 their 
Fig. 9) investigated the relationship between Li abundances and masses for 
ZAMS stars in the Pleiades clusters, and found that, when the mass of the 
stars ranges between 0.9M o to 1.4M , the Li abundances do not strongly de- 
pend on the mass. We selected a sample of WTTS stars in the range of mass 
0.9-1.4M Q in Taurus-Auriga star forming regions (SFRs, in a same nebular, 
star has the similar metallicity). 



We concentrate on the connection between the EW(Li) and P ro t- First, since 
the Li I A6707 line is very strong in T Tauri stars, it is on the saturated part of 
the curve of growth, so small errors in the mea s ured Li I strength will lead to 
large error in the derived abundances ((Duncanl . Il991l ). Using just the EW(Li) 
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can avoid some uncertainties propagated from uncertainties of T e ff estimates 
which are provided by the estimates from (B-V) and (V-I) colors. Second, 
the measurements of rotation velocity v sin i yields only a lower limit to the 
rotation velocity due to the unknown angle of inclination, i. We use both 
the directly observed EW(Li) and P rot instead of Li abundance and rotation 
velocity of stars in our rot at ion- lithium study. 

This paper is organized as follows. The periods and Li data of WTTSs are 
presented in Sect. 2. The relationship between EW(Li) and P ro t and discussion 
provided in Sect. 3. And the results are summarized in Sect. 4. 



A The Rotation periods and Li data 



We selected from the literature a sample of WTTSs which according to their 
position in the H-R diagram are in the mass range from O.9M to 1.4M , 
and whose EW(Li) and P rot have been determined. The location of the stars 
in the H-R diagram is shown together with theoretical pre-main sequence 
evolutionary tracks in Fig. 1. The theo r etical pre-main sequence evolutionary 
tracks were taken from Cohen fc Kuhi ( 1979f ). The known binaries have been 
rejected, since the late-type binary system may be tidall y locked rotation 



which lead naturally to slower lithium destruction rates ((Maccarone et al 
2005b . 



A.l Rotation periods 



Our sample consists of 21 WTTSs in the range of mass O.9-1.4M in Taurus- 
Auriga. In table 1, the up-to-date rotation periods and equivalent widths are 
presented for our sample stars. The rotation periods of these stars are taken 
directly from the following literature: 



The periods of WTTSs from Bouvier et al. ( 19931 ) are shown table 1. Our 
sample includes the objects with a multi-site pho tometric campaign to mon- 



itor T Tauri stars over more than two month bv lBouvier et al. (1993). Two 



methods, namely the string-length method and periodogram analysis, were 
used to derive the P rot for ten WTTS stars, mo st of them with confi dence 
level larger than 99%. The secondary sample of Bouvier et al. ( 19931 their 
Table 4) is also included. Totally, we selected seven WTTSs from their sam- 
ple. 

Xing et aD ( 2006f ) selected a sample of X-ray sources that have been identi- 



fied as WTTS stars around the Taurus- Auriga SFRs. They monitored the 
light variations for 22 WTTSs and obtained the P rot for 12 stars using two 
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methods: the Phase Dispersion Minimization (PDM) and Fourier analysis 
methods (PERIOD04), most of them with confidence level large than 99%. 
We took seven WTT Ss from their sample. 

Bouvier et all f|l997h monitored the light variations of 58 WTTSs and de- 



rived photometric periods for 18 stars using 3 methods: the periodogram 
analysis, CLEAN deconvolution algorithm and string-length estimator. Ex- 
cept for RXJ0409. 2+2901, all periods are detected at the 99% confidence 
level in the per iodogram. Five objects have been taken from their sample. 
Grankin (1993) presented the results of BVR photometry of 22 WTTSs. 



They obtained rotation periods for ten stars. Two WTTSs were selected 
from this sample. 



A . 2 Equivalent widths of Li 



The equivalent widths of Li for the above stars are taken directly from the 
following literature: 



The EWfLi) of the iBouvier et al.l (|l993h sample were taken directly from 



Basriet et al.l (|199lh . They reported the observations of strong Li I A6707 



line in 28 T Tauri stars in the Taurus- Auriga star formation complex. Line 
strengths were obtained using high resolution spectra from the Hamilten 
echelle at Lick Observatory. They have corrected the Li I equivalent widths 
for continuum veiling based on a simultaneous measurement of the actual 
veiling present. 

The EW(Li) of the lBouvier et all (|l997l ) sample a nd three stars fRX J0430.8- 
RX J0405. 1+2632 and RX J0432.7+1853 ) of th e lXing et al.1 (|2006l ) sample 
were taken directly from lWichmann et al.l (J2000J). They presented a detailed 
study for Li-rich stars discovered by ROSAT in Taurus-Auriga SFRs 
results are based on high-resolution echelle sp e ctra. 



-2113, 



the 



spe ctra. 

The Li I equivalent widt hs of the Granki 3 (|l993h 



NTTS045251+3016 from IXing et alT l2006) 



were 



Walter et al. ( 19881). measur e d from high-dispersion spectra. 



sample and the star 
taken from the list of 



Two stars of the IXing et all fcOOfil) samp le, [L#98]37 and [LF98]53, were 
identified as WTTSs bv lLi fc HuT ((19981 ) . and the Li I equivalent widths 
were taken from their results, which are based on t he intermediat e resolution 
spectra. The EW(Li) of HD 287927 was taken from lWalterl (|l986l ). also based 
on the intermediate resolution spectrum. 
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B Results and discussion 



B. 1 Li-rotation for WTTS 



We plotted the EW(Li) versus P rot for our sample stars in Fig. 2. This shows 
that there is a clear correlation between the EW(Li) and P rot , i.e. the rapid 
rotators have lower EW(Li) and the depletion of lithium proceeds fast, when 
the rotation period of star evolves towards 3 days for WTTSs with mass 
between 0.9M Q and 1.4M in Taurus-Auriga SFRs. 



Our results are in good agreement with the hypothesis that surface Li depletion 
takes place during PMS evolution for low-mass stars as a r esult of Li burnin g 
via (p, a) reactions at low temperature s of T > 2.6 x 10 6 K dKing et aL 2.000). 
These results support the conclusion by Piau fc Turck-Chieze I ([20021 their Fig. 
4) about a model for lithium depletion with age of star. Their results predict 
that the surface Li depletion will take place during PMS evolution for low-mass 
stars. They also distinguish two phases in lithium depletion: (1) a rapid nuclear 
destruction in the T Tauri phase before 20 Myr whatever the mass in the range 
between 0.8 and 1.4M , and (2) a second phase where the destruction is slow 
and moderate, which is largely dependent on the hydrodynamic instability 
located at the base of the convective zone. 



It is important to know whether the EW(Li) depends on the masses or not. 
We plotted the EW(Li) against the masses for our sample stars in Fig. 3. This 
indicates that the EW(Li) of our sample stars does not show any trend as a 
function of mass. It means that the EW(Li) do not depend on the mass of 
WTTSs for our sample stars. 



In order to further confirm our result, we plotted the EW(Li) against P ro t 
for WTTSs i n Orion in F i g. 4. T he EW(Li) and P rn t of these stars are taken 
directly from Alcala et al. ( 1996f ) and Marilli et al. ( 20051 ) . These selected stars 
are late G or early K-type stars whose masses are close to the masses of WTTS 
sample in Taurus-Auriga (since the mass is not direct observable, except in 
eclipsing b inaries, i t wou ld be better to use T e ff or spectral type to estimate 
the mass (|Martml . Il997h . In the absence of a luminosity determination for 
these stars, the range of mass was estimated from the spectral type of the 
stars). Fig. 4 shows that the correlation between EW(Li) and P rot of WTTSs 
(late G and early K-type) in Orion nebula is in a good agreement with that 
one of WTTSs (O.9M < M < 1.4M ) in Taurus-Auriga. It is evident that the 
rapidly rotating stars have lower EW(Li), although the stars in Taurus- Auriga 
differ from that in Orion in some aspects. 
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B.2 Li-rotation for young solar type stars 



The l ithium-rotation relation for ZAMS stars in Plei ades (e.g.|_T schape & Riidiger 
2001) and young low- mass man sequence stars (e.g. Rebolo &: BeckmanT 19881 : 
Chaboverl . Il998f ) has long been known. In "older" clusters (Pleiades, Hyades ) 
the faster rotating stars show less the Li depletion (jTschape fc RudigerLl200lh . 
In order to compare the possible effect of rotation upon lithium depletion 
between young solar- type mai n sequence stars and PMS, we take the P rn t . 
and EW(Li) for 7 ZA MS stars ( Krishnamurthi et allll9980l : ISoderblom et al. 



1993; iMessinal. l200lh in the Pleiades cl uster and 12 young low-mass main se 
quence stars ( Rebolo fc Beckmanl . 19881 ) in the Hyades cluster, and plotted the 
EW(Li) versus P rot of these stars with WTTSs in our sample in Fig. 5. The 
spectral types of these young stars were in the range between late G and early 
K-type (in the r ange of mass 0.9- 1 .4Mq , determined from mass-temperature 
relation given bv lSoderblom et al.l l)l99ffl ). The EW(Li), P rot and other stellar 
properties of these stars in Pleiades and Hyades clusters are shown in table 3. 

Fig. 5 shows that the turn over of the lithium-rotation relation curve is at 
ZAMS phase. The EW(Li) decreases with decreasing P rot when stars young 
than ZAMS, whereas EW(Li) decrease with increasing P rot (decreasing rota- 
tional velocity) when stars are on the ZAMS or older. 



These results are in good agreement with the rotation evolution model (e.g. 
Bouvier et al. . 1994 : Soderblom et al. . 19931 : Cemeron et al. . 19951 : Keppens et al 
'., 995) a nd Li could serve as a "clock" of stellar evolution in the PMS phase 
( Drake! liiolh . The rotation evolution model and the results of these obser- 
vations indicate that the main features of the rotation evolution of low-mass, 
late type stars are the strong PMS spin up from moderate rotation in the T 
Tauri phase to ultrafast rotation at ZAMS and an increase of the rotation 
period of star (spin down) with increasing age for main sequence stars. 



The relation between EW(Li) and P rot in Fig. 5 also shows that rapid rotators 
have lower EW(Li). This result is consistent with the lithium-age correlation 
in the sense that there is less lithium in the surfaces of older WTTSs (Post- 
TTSs) than in young WTTSs. e.g. the EW(Li) of older WTTSs (TAP 9 whose 
location in H-R diagram very near ZAMS star) is lower than that of young 
WTTSs (TAP 57 whose location in H-R diagram still on Hayashi line). 



C Summary 



In this work we have discussed the correlation of lithium-rotation. Our main 
conclusions can be summarized as: 
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• At least for WTTSs with mass between O.9M and 1.4M in Taurus-Auriga 
Nebula, there is a clear correlation between EW(Li) and P ro t, i.e. on aver- 
age, rapidly rotating stars have lower equivalent widths of Li. That can be 
explained as an age effect of Li depletion during pre-main sequence. It is 
clear that rapid rotation does not inhibit Li depletion among low mass PMS 
stars. 

• The fast burnings of Li begin at the phase when the rotation period of the 
star evolves to approach 3 days. And the surface lithium depletion always 
happens during the PMS phase. 

• The turn over of the lithium-rotation connection curve at the phase of 
ZAMS. The equivalent widths of Li decreases with decreasing rotation pe- 
riod when stars are younger than ZAMS, whereas the equivalent widths of 
Li decreases with increasing rotation period (decreasing rotation velocity) 
when stars evolve beyond the ZAMS. 
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Table .1 

Rotation al periods a nd st ellar propertie s for 21 WTT Ss. Referen c e to table 



pr 

1: XfcL: IXinget all EM) and Hi J Jlflflft: XfcW: king et all jB) and 

Wichma nn et al.| (1200(1): XfcWa: IXing et alJ \2Q0(j) an d IWalterl dl98fih: XfcG : 



G&W: 



Xing et ahlfcoOpl i andlGregorio-Hetermet Heteml fcOQ^ : R:lRonvier Rt allilflfl^: 
" Grankinl (Il996h and IWalter et all (|l988j ): R&W: iRouvier et all <ll997h and 



Wichmann et al.l (|2000l ). The mass of stars are taken directly from ab ove literature 



or 



from comparison with evolutionary tracks of lCohen Kuhl {l979). 



Star 


Pr oi (d) 


log(L*/L ) 


SpT. 


T e// 


EW(Li)(mA) 


M(M ) 


ref 


The sample of Xing et al. photometry 


[LH98}37 


1.13 


0.12 


K0IV 


5236 


240 


1.2 


X&L 


[LH98}53 


0.728 


-0.05 


G2IV 


5792 


230 


1.0 


X&L 


HD 287927 


0.772 


-0.12 


G5 


5554 


200 


0.98 


X&G 


NTTS 045251+3016 


9.12 


0.04 


K5 


4034 


580 


0.9 


X&Wa 


RX J0405. 1+2632 


1.93 


-0.34 


K2 


4897 


219 


0.94 


X&W 


RX J0430.8+2113 


0.741 


0.19 


G8 


5309 


141 


1.27 


X&W 


RX J0432.7+1853 


1.55 


-0.1 


Kl 


5105 


253 


1.1 


X&W 


The sample of Bouvier et al. (1993) photometry 


V1068 Tau 


3.37 


0.04 


K7 


4060 


510 


0.9 


B 


V836 Tau 


7.0 


-0.22 


K7 


4060 


570 


0.9 


B 


NTTS 045226+3013 


2.24 


0.17 


K0 


5240 


440 


1.2 


B 


NTTS 034903+2431 


1.6 


-0.3 


K5 


4395 


370 


1.03 


B 


NTTS 041636+2743 


5.64 


-0.04 


K7 


4060 


600 


0.9 


B 


RX J043005. 1+181351 


2.7 


0.11 


K2 


4950 


420 


1.24 


B 


RX J043214.9+182013 


3.75 


0.045 


K7 


4060 


570 


0.9 


B 



The sample of Bouvier et al. (1997) photometry 



RX J0409. 2+2901 


2.74 


0.0 


Kl 


5105 


413 


1.2 


B&W 


RX J0415. 4+2044 


1.83 


0.0 


K0 


5236 


270 


1.1 


B&W 


RX J0423.7+1537 


1.605 


-0.2 


K2 


4855 


361 


1.0 


B&W 


RX J0438.7+1546 


3.07 


0.1 


Kl 


5105 


419 


1.2 


B&W 


RX J0457.2+1524 


2.39 


0.2 


Kl 


5105 


446 


1.4 


B&W 






The sample 


of Grankin (1993) 


photometry 






NTTS 041559+1716 


2.52 


-0.4 


K7 


4438 


530 


0.93 


G&W 


NTTS 042835+1700 


1.55 


-0.52 


K5 


4352 


150 


0.88 


G&W 
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Table .2 

Rotation periods and stellar properties for 9 s tars in P l eiade s clu sters and 14 stars 



in Hya des clusters. Reference to table 2: MfeS: lMessinal (|200lJ ) and lSoderblom et al 
(1993); R: iRebolo fc Beckmanl (|l988h 



Pleiades clusters stars 



HII 




SpT. 


T eff 


EW(Li) 


ref 


number 


(d) 




°K 


mA 




263 


4.82 


G8V 


5060 


290 


M&S 


345 


0.84 


G8V 


5160 


245 


M&S 


738 


0.83 


G9V 


5140 


203 


M&S 


1039 


0.784 


K2V 


4720 


333 


M&S 


2244 


0.56 


K2.5V 


4720 


268 


M&S 


882 


0.581 


K3V 


4500 


212 


M&S 


1883 


0.235 


K2V 


4560 


282 


M&S 


3197 


0.44 


K3v 


4440 


302 


M&S 


1653 


0.74 


K4.5Ve 


4220 


108 


M&S 


Hyades clusters stars 


BD+ 


Perioc 


SpT. 


T eff 


EW(Li) 


ref 




(d) 




°K 


mA 




19 694 


9.2 


G5 


5460 


32 


R 


17 707 


9.1 


G8V 


5570 


29 


R 


18 623 


5.5 


GOV 


6060 


86 


R 


16 589 


6.2 


GOV 


5930 


72 


R 


16 592 


7.9 


G2V 


5840 


69 


R 


18 636 


6.1 


G5V 


5650 


70 


R 


15 624 


5.1 


GO 


6070 


84 


R 


16 601 


8.5 


G2V 


5770 


51 


R 


15 627 


5.9 


GO 


6200 


85 


R 


16 606 


7.4 


GV 


5920 


82 


R 


17 731 


3.2 


GO 


6340 


19 


R 


17 734 


11.4 


G5 


5230 


3 


R 


15 642 


9.0 


G5 


5540 


15 


R 


15 651 


6.5 


GO 


5940 


84 


R 
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Period (days) 

Fig. .2. Lithium A6707 equivalent widths of WTTSs (mass range 
O.9M0 < M < IAMq in Taurus-Auriga SFRs) as a function of the rotation 
periods. 
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CO 



0.9 



1.1 
M(M ) 



1.2 



Fig. .3. Lithium A6707 equivalent widths of WTTS (mass range cover by 
O.9M < M < IAMq in Taurus-Auriga SFRs) as a function of their masses. 
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Fig. .4. Lithium A6707 equivalent widths for WTTSs in Orion nebula as a function 
of their rotational periods. The spectral type of sample stars ranges from G9 to K5. 
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Fig. .5. The lithium A6707 equivalent widths as a function of the rotational periods 
of WTTSs of mass 0.9M o < M < 1.4M in Taurus-Auriga SFRs (crosses), ZAMS 
stars between late G and early K-type in Pleiades cluster (filled square) and young 
solar-like stars between late G and early K-type in Hyades cluster (stars). 
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